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INTRODUCTION 

This program investigates the fluid dynamics and 
simultaneous solidification of molten solder droplets 
impacting on a flat substrate. The problem of interest 
is directly relevant to the printing of microscopic sol- 
der droplets in surface mounting of microelectronic 
devices. The study consists of a theoretical and an 
experimental component. The theoretical work uses 
axisymmetric Navier-Stokes models based on finite 
element techniques. The experimental work will be 
ultimately performed in microgravity in order to allow 
for the use of larger solder droplets which make feasi- 
ble the performance of accurate measurements, while 
maintaining similitude of the relevant fluid dynamics 
groups (Re, We). 

A schematic of the two distinct stages (flight and 
impact) of the problem examined is shown in Fig. 1. 
The primary application of interest (solder microdro- 
plet dispensing) employs solder droplets approxi- 
mately 50 to 100 pm in diameter, which collide, 
spread, recoil and eventually solidify on the substrate. 
Due to the small size of the droplets and the relatively 
high surface tension coefficient of solder, gravity ef- 
fects are negligible. This solder application technol- 
ogy has shown great promise in microelectronic pack- 
aging and assembly, therefore, the development of a 
good understanding of the pertinent fluid dynamics 
and solidification phenomena is essential for its suc- 
cessful commercial implementation. However, prog- 
ress in this area has been hindered by the small length 
scales of the problem (50 to 100 pm), which have 
made experimental measurements of the relevant 
transport phenomena difficult. Alternative ap- 
proaches, which employed much larger (mm-size) 
droplets, yielded results that were affected by the 
masking effects of gravity. Hence, even though mm- 
size droplets yield significantly improved resolution, 
the applicability of the obtained results for much 
smaller droplets remains suspect. Conducting ex- 
periments in a microgravity environment eliminates 
the unwanted influence of gravity and makes the ex- 


perimental investigation of large droplet dispension 
directly relevant. 


The Reynolds, Weber and Froude numbers char- 
acteristic of the process shown in Fig. 1 are defined by 
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where V 0t d 0 denote droplet impact velocity and di- 
ameter, while p, v, /correspond to the density, kine- 
matic viscosity and surface tension coefficient of the 
liquid. To exemplify the disparity in the importance of 
gravity in the dispension of large and small droplets, 
the values of Re, We and Fr were calculated for a set 
of parameters corresponding to the real dispension 
process in normal gravity (50 micron diameter solder 
droplet impacting on a flat surface with a velocity 
lm/s). These values were Re=157, We=1.2, Fr=2038, 
and illustrate the importance of inertia and surface 
tension, as well as the insignificance of gravity effects. 
If a larger (1mm) drop is used in normal gravity ex- 
periments, a slower impact velocity V 0 is required for 
similitude based on Re. In turn, similitude in terms of 
the Weber number requires a smaller surface tension 
coefficient. To maintain similitude in this specific 
example, the surface tension coefficient as well as the 
impact velocity of the larger droplets need be reduced 
by twenty-fold. The resulting values of the above di- 
mensionless groups are then Re=157, We=1.2 and 
Fr=0.25. Clearly, the drastic decrease in the value of 
the Froude number proves that gravity effects become 
significant for mm-size droplets, and that the presence 
of a microgravity environment is necessary in large 
solder-droplet impact experiments. 


OBJECTIVE 

The study aims to create a science base and iden- 
tify the influence of the dominant process parameters 
in solder droplet dispensing. These parameters are: 
droplet size and velocity; droplet, substrate and ambi- 
ent gas temperatures; and contact angle between sol- 
der and substrate before and after solidification. The 
sensitivity of the solidified-droplet (bump) shape and 
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size to variations in the above parameters is critical 
because solder bump volume, position, and height 
variation are key metrics for solder jet technology. 
Through a combination of experiments and numerical 
modeling, the effect of the dimensionless groups de- 
fined in Eq. (1) and the physics they represent are 
systematically documented. 

METHODOLOGY 

The research consists of a theoretical and an ex- 
perimental component. The theoretical component 
investigates the fluid dynamics and solidification of a 
molten solder droplet during its impact on the sub- 
strate, in order to attain a fundamental understanding 
of the miniature solder deposition process. The ex- 
perimental component tests the numerical predictions 
and provides necessary input data (such as wetting 
angles) for the theoretical model. Details of both com- 
ponents are given below. 

Theoretical 

The model simulates the axisymmetric impact and 
subsequent solidification of an initially spherical, 
molten solder droplet on a flat multi-layer composite 
substrate. When a gravitational field is present, the 
droplet is injected along that direction, which is or- 
thogonal to the target surface. The Navier- Stokes 
equations combined with heat transfer and solidifica- 
tion are solved in the liquid phase using a Lagrangian 
approach. The heat conduction equation is solved in 
the solid phase, i.e., substrate and solidified sector of 
the droplet. It is noted that the employed axisymmetric 
approach is valid for situations where the solder 
droplet impinges on the substrate at an angle equal or 
approaching 90°. The mathematical model formula- 
tion has been described elsewhere in detail (see Wald- 
vogel et al., 1996; Waldvogel and Poulikakos, 1997) 
and will not be repeated here for brevity. However, a 
brief description is given in the following for com- 
pleteness. 

Fluid Dynamics: Laminar flow of a constant property 
fluid (molten solder) is assumed. The radial and axial 
components of the momentum equation are considered 
along with the continuity equation. The artificial 
compressibility method is employed, thus introducing 
a pressure term to the mass conservation equation. 
The free surface deformation is tracked by the two 
corresponding stress balance equations (along r and 
z). The wetting force at the dynamic contact line be- 
tween the liquid droplet and the substrate is neglected 


throughout the analysis. This assumption is valid 
during the initial stages of droplet impact where iner- 
tia forces are high with respect to the retarding forces 
of viscosity and surface tension. In the case of solder 
jetting, solidification occurs very quickly at the contact 
area, thus eliminating the need for a wetting condition 
altogether. The moving contact line is treated by al- 
lowing for slipping along the radial direction (Dussan, 
1979). 

Heat Transfer: The Lagrangian formulation is used 
for the conjugate heat transfer process in the droplet 
and the substrate. The energy conservation equation is 
solved in the droplet interior as well as in the sub- 
strate. Heat transfer from all exposed surfaces is ne- 
glected. 

Thermal Contact Resistance: Thermal contact re- 
sistance between droplet and substrate is modeled by a 
thin layer of arbitrary thickness. This layer is assigned 
zero heat capacity and experiences only axial conduc- 
tion. Thermal contact resistance is quantified by a 
dimensionless coefficient (Biot number), defined by 


where h c denotes the heat transfer coefficient for the 
thin contact area between the droplet and the sub- 
strate, k 0 the thermal conductivity of the solder and d 0 
the droplet diameter before impact. Unfortunately, no 
correlations published to date have been identified to 
determine the value of h c in the parametric domain of 
solder jetting. To this end, specific attention is given 
to identify the thermal contact resistance in solder 
microdroplet dispensing and, in turn, facilitate future 
studies in this area. 

Solidification: The following assumptions are made 
regarding the solidification of the molten-solder 
droplet. First, since a eutectic alloy is studied, a sharp 
boundary separating distinct liquid and solid regions 
is included in the model. Furthermore, it is assumed 
that phase change occurs at the equilibrium freezing 
temperature (i.e., no undercooling is considered). It is 
noted that the droplet solid and liquid phase densities 
are nearly identical since the 63% Sn-37% Pb solder 
material experiences only a 0.2% contraction upon 
freezing. The exact specific heat method proposed by 
Bushko and Grosse (1991) was adopted to model so- 
lidification. This approach facilitates the finite ele- 
ment formulation and has the added advantage that it 
is very accurate in conserving energy as the droplet 
solidifies. 
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Numerical Solution Procedure 

The computational domain is discretized with a 
mesh of triangular elements, and the numerical model 
is solved using a Galerkin finite element method. 
Quadratic and linear shape functions are used for ve- 
locity and pressure, respectively. An implicit method 
is utilized for the numerical integration of the fluid 
dynamics equations in time, while a Crank-Nicholson 
scheme is used for the energy equation. The details of 
the iterative solution technique and mesh generation 
algorithm are given in Waldvogel and Poulikakos 
(1997). The spatial resolution needed to provide grid 
and time-step independence is '■600 nodes in the 
droplet and ~ 1 , 1 00 nodes in the substrate. A dimen- 
sionless time step of 5t V„ /d 0 = 5xl<T 4 was found ade- 
quate for time-step insensitive results. 

Experimental 

Two separate jetting apparati are employed in the 
experiments. The first is capable of jetting 50-80|im 
solder droplets with velocities ~lm/s, while the other 
generates mm-sized droplets at injection speeds of 
~0.1m/s. Both devices create droplets on demand, and 
feature heated reservoirs, in which high-purity 
63%Sn-37%Pb solder is maintained at temperatures 
~210°C (melting point of solder is 183° C). The first 
device uses nitrogen gas to pressurize the system 
forming the droplets, thus hindering the formation of 
oxides within the reservoir. Individual droplets are 
generated as a voltage pulse excites a piezoelectric 
transducer mounted within the device. The second 
device (which forms larger droplets) utilizes mechani- 
cal forcing of the solder pool (via a plunger) through a 
round orifice at the bottom of the reservoir. The drop- 
let flight occurs in a nitrogen atmosphere to avoid 
oxygen adsorption on the solder surface. The presence 
of oxygen is avoided because it is known to degrade 
the surface properties of liquid solder. The impact 
velocity V 0 of the small droplets is determined using a 
stroboscope and a steady stream of droplets generated 
at a known injection frequency. For the larger drop- 
lets, a high-speed video system is used to measure 
droplet impact velocities. 

The solidified shape of the small deposited droplets 
is evaluated using scanning electron microscopy 
(SEM). The experimental results reported in the fol- 
lowing section were all obtained in normal gravity. 
Small droplets impacted on a two-layer composite 
substrate (Ni on Si). For the larger (mm) droplets, the 
substrates consisted of a single-layer metal material. 


RESULTS AND DISCUSSION 

Thermal Contact Resistance of Microdroplets 

The thermal contact resistance between a typical 
electronic substrate and impacting solder droplets was 
examined first. The chosen conditions correspond to a 
63%Sn-37%Pb solder droplet with pre-impact di- 
ameter </a=53(J.m, velocity Vq= 1.6 m/s and temperature 
2 1 0°C. The two-layered substrate (consisting of a 
2|im-thick nickel layer on a 212pm layer of silicon) 
had an initial temperature of 35°C. The dimensionless 
numbers corresponding to these baseline conditions 
are Re=266, We=3.23 and Fr=4930. These operating 
conditions are characteristic of practical situations, 
and have been chosen primarily because of the avail- 
ability of experimental data; see Waldvogel and Pou- 
likakos (1997), as well as Xiong et al. (1998). 

A series of model simulations was performed with 
all process parameters fixed to the baseline values 
stated above, while the values of Bi were varied over a 
range from 0.01 to 0.75. These values of Bi corre- 
spond to contact heat transfer coefficients h c in the 
range 4-350 kW/m 2 K. Even though the value of Bi 
changed from run to run, it was kept constant 
throughout each of these simulations. The computed 
solder bump outline for each value of Bi revealed that 
this shape is very sensitive to the contact resistance 
between droplet and substrate; see Xiong et al. (1998). 
Figure 2 presents two geometric parameters describing 
the solidified-bump shape and its dependence on 
thermal contact resistance. Both parameters are plot- 
ted in their dimensionless form versus Biot number 
(likewise, heat transfer coefficient h c ). Two values of 
height (z/do) for each solidified bump appear in Fig. 2: 
one on the symmetry axis and another at the off-axis 
location of maximum elevation. It is clear that for all 
values of Bi examined, these two heights are essen- 
tially identical. Figure 2 also shows that the predicted 
bump height varies by up to 20% depending on ther- 
mal contact resistance; the maximum height is real- 
ized at Bi= 0.1 (or *,=47.2 kW/m 2 K). These variations 
of bump dimensions with thermal contact resistance 
can be critical in solder jetting, as solder bump height 
deviations of only a few percent can be detrimental to 
the successful implementation of the technique in 
high-precision electronic component manufacturing. 
The dimensionless contact area ( d c0 Jd 0 f between the 
solidified solder droplet and the substrate is also 
shown in Fig. 2 and demonstrates a weaker depend- 
ence on thermal contact resistance, with an overall 
increasing trend towards higher values of Bi or h c . 
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According to this trend, more efficient heat transport 
results in larger contact area, and, in turn, better at- 
tachment of the bump on the substrate. 

A detailed comparison of several computed bump 
shapes with those determined experimentally was 
conducted, see Xiong et al. (1998), to determine the 
value(s) of Bi which provided the best match between 
the two shapes. Shape comparisons were conducted by 
means of the morphology of the ripples on the outline 
of the bump (see Fig. 3), and by measuring the contact 
area at the bump base. A phase-dependent value of Bi 
allowed best agreement between experiment and mod- 
eling. More specifically, the value Bif= 0.5 ( h c - 235 
kW/m 2 K) was identified for liquid solder/substrate 
contact, and the value Bi s = 0.0275 ( h c - 13 kW/m 2 K) 
for solid solder/substrate contact. Figure 3 a shows the 
bump geometry obtained in the experiments under the 
baseline conditions. Figure 3b displays the respective 
model-predicted shape when the above-mentioned 
values of Bi\ and Bi s were used. It was concluded that 
the best agreement between experiment and modeling 
occurs when the thermal contact between solder and 
the substrate is reduced nearly by twenty-fold with the 
appearance of the solid phase. 

Impact of Large Droplets 

The results presented in this section were obtained 
in normal gravity and are discussed to demonstrate the 
capabilities of the solder-droplet generation apparatus 
and the high-speed visualization equipment. Both will 
be implemented in the pending microgravity experi- 
ments which will involve mm-sized droplets. 

A 63%Sn-37%Pb solder droplet with pre-impact 
diameter ^=0.9inm, velocity Fo=0.6m/s and tem- 
perature 210°C was visualized using high-speed video 
during its collision with a smooth stainless steel sub- 
strate at a temperature of 20°C. The dimensionless 
numbers corresponding to these conditions are 
Re=1676, We=7.4 and Fr=37. Figure 4 presents the 
temporal variation of two geometric parameters de- 
scribing instantaneous droplet shape from impact up 
to complete solidification. These parameters are 
maximum height z max and maximum splat diameter 
d maXi both of which change during the dynamic stages 
of the impact event. Both maximum height and di- 
ameter have been normalized in Fig. 4 with respect to 
d 0 \ time is also presented in a dimensionless form, i.e., 
t*=tVo/d 0 . The oscillatory nature of the event is appar- 
ent in both curves of Fig. 4. The error bars drawn over 
selected points of the curves shown in Fig. 4a and 4b 
indicate the uncertainty tied to pixel size resolution. It 


is noted that time zero in Fig. 4 corresponds to the 
first acquired image where contact of the droplet and 
substrate was visualized. In this case, time zero does 
not correspond to the instant when contact between 
the two bodies was first established (time resolution is 
lms). The shape of the solidified bump is shown by 
the inset of Fig. 4b. 

Figure 5 presents selected experimental images 
obtained via high-speed video when a 0.9mm solder 
droplet impacted on a smooth copper substrate of tem- 
perature 25°C. The initial droplet temperature was 
210°C, and the impact velocity was 1.14m/s. The non- 
dimensional numbers corresponding to these condi- 
tions are Re - 3200, We — 27.7 and Fr = 147.5. In the 
shown sequence, time proceeds from left to right and 
from top to bottom. Under the above conditions, 
breakup of the droplet was revealed at the end of the 
first recoiling cycle after impact. The smaller droplet 
created by this breakup was subsequently reattached to 
the main mass, eventually creating a solid bump 
showing no apparent signs of the earlier separation 
event (see bottom right frame). 
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Fig. 1 Basic stages in solder microdroplet dispensing: (a) Flight, (b) Impact. 
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Fig. 2 Geometric parameters describing the solder bump shape and its dependence on thermal contact resistance: 
z/do denotes the dimensionless height either on the symmetry axis (circles) or at the off-axis location of maximum 
elevation (squares); (dc 0n /do) 2 denotes the dimensionless contact area between the solidified solder droplet and the 
substrate. 



Fig. 3 (a) Solder bump geometry obtained in the microdroplet experiments for the baseline jetting conditions 
(Re=266, We=3.23 and Fr=4930). (b) Model -predicted shape providing best agreement between experiment and 
modeling. In order to achieve this agreement, a phase-dependent value of Bi was used: 5i/=0.5 for liquid sol- 
der/substrate contact, =0.0275 for solid solder/substrate contact. 
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Fig. 4 Temporal variation of (a) dimensionless maximum droplet height z max /do, and (b) maximum splat diameter 
d m ax/do from impact up to complete solidification of a mm-sized solder droplet impacting on a stainless steel sub- 
strate (Re=1676, We=7.4 and Fr=37). The error bars indicate the uncertainty due to pixel size resolution. The 
shape of the solidified bump is shown in the inset of frame (b). 



Fig. 5 Selected experimental images obtained using high-speed video when a 0.9mm solder droplet impacted on a 
smooth copper substrate (Re = 3200, We = 27.7 and Fr = 147.5). Time proceeds from left to right and from top to 
bottom. Under these conditions, breakup of the droplet was observed at the end of the first recoiling cycle after 
impact. 


239 




